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Large mass bolometers are used in particle physics experiments to search for rare processes. By operating 
at low temperature, they are able to detect particle energies from few keV up to several MeV, measuring 
the temperature rise produced by the energy released. This study was performed on the bolometers of the 
CUORE experiment. The response function of these detectors is not linear in the energy range of interest, and 
it changes with the operating temperature. The non-linearity is found to be dominated by the thermistor and 
its biasing circuit, and is modeled using few measurable parameters. A method to obtain a linear response is 
the result of this work. It allows a great simplification of the data analysis. 
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I. INTRODUCTION 

Bolometers are detectors where the energy of parti- 
cle interactions is converted into heat and measured via 
temperature variation. They provide an excellent energy 
resolution, although they are slow compared to conven- 
tional detectors. These features make them a suitable 
choice for experiments looking to rare processes, like neu- 
trinoless double beta decay and dark matter interactions. 

The CUORE experiment searches for neutrinoless dou- 
ble beta decay in 130 Tei, using bolometers made of 
Te02 crystals whose temperature is measured by neutron 
transmutation doped (NTD) Germanium semiconduc- 
tor thermistors-Sal. Operated at temperatures of about 
lOmK, these detectors maintain a resolution of few keV 
over their energy range, extending from few keV up to 
several MeV. In this range the response function is found 
to be non-linear. The conversion from signal amplitude 
to energy is complicated and the shape of the signal de- 
pends on the energy itself. Moreover the amplitude of 
the signal depends on the temperature of the detector, 
which cannot be kept stable by current cryostats at a 
level that would not perturb the resolution. Understand- 
ing the non-linearities is important to have a reliable cal- 
ibration, to control the signal shape and to improve the 
resolution. 

Electrothermal models for large mass bolometers have 
been already proposed^—. Based on systems of differen- 
tial equations, they describe the behavior of the detector 
and reproduce the shape of the signal and some features 
of the calibration function. However, they rely on many 
parameters that are difficult to measure with the neces- 
sary precision in working conditions. These models are 
useful to describe the response of a bolometer, but they 
cannot produce tools to be used on a real detector. 

In this paper a model based on few parameters is pro- 
posed. It shows that the non-linearities are dominated 



by the thermal sensor and by the biasing circuit that is 
used to read it. These elements are fully measurable in 
working conditions, so that the model can be applied to 
the detector. A method to remove the non-linearities is 
derived. Its application in the data analysis linearizes the 
response function, improving the results. 



II. DESCRIPTION OF THE DETECTOR 

A CUORE bolometer is composed of two main parts, 
a T3O2 crystal and a NTD-Ge thermistor. The crystal is 
cube-shaped, weighs 750 g and its heat capacitance C at 
the working temperature of 10 mK is of order 2nJ/K 7 . 
It is held by Teflon supports on copper frames, which 
are connected to the mixing chamber of a dilution refrig- 
erator. The thermistor is glued to the crystal and the 
biasing wires are glued to the copper frames (see Fig. [T]) 
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Figure 1. Sketch of a CUORE bolometer. The Te0 2 crystal 
is held by Teflon supports, the thermistor is glued to the crys- 
tal and its wires are glued to copper frames. The supports 
and the thermistor wires thermally couple the crystal to the 
copper frames, which act as heat bath. 
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When an amount of energy E is released in the crystal, 
its temperature increases by E/C and then returns to 
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equilibrium via the thermal conductance of the supports 
and of the thermistor. The latter acts as thermometer, 
converting temperature, T, into resistance, R, according 
to the relationship^: 



R(T) = R exp (To/T) 1 



(1) 



where i?o and To are parameters that depend on the di- 
mensions and on the material of the thermistor. Their 
value is about 1.1 S7 and 3.4 K respectively. At lOmK 
the parameter 7 can be considered constant and equal to 

1/2 aia. 

To read out the signal, the thermistor is biased with 
constant current, which is provided by a voltage genera- 
tor, Vb, and a load resistor, R L , in series with the ther- 
mistor. The resistance of the thermistor varies in time 
with the temperature, R(t), and the voltage across it, 
V(t), is the bolometer signal. The value of Rl is cho- 
sen to be much higher than R(t) so that V(t) is almost 
proportional to R(t). Since from Eq. (fTJ) positive tem- 
perature variations induce negative resistance variations, 
the voltage Vb is chosen to be negative in order to obtain 
positive signals. The connecting wires add in parallel to 
the thermistor a parasitic capacitance c p (see Fig. 
The signal V(t) is amplified, filtered with a 6-pole active 
Bessel filter and then acquired with an 18 bit DAQ sys- 
tem. The front end electronics, which provide the bias, 
the load resistors, and the amplifier, are placed outside 
the cryostat, at room temperature^. 
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Figure 2. Biasing circuit of the thermistor. A voltage gener- 
ator Vb biases the thermistor resistance R(t) in series with a 
load resistance Rl- The bolometer signal is the voltage V(t) 
across R(t). The wires used to read V(t) have a non- negligible 
capacitance c p . 



At 10 mK the value of R(T) is of order 100 M O, R L 
is chosen as ~ 50Gf2 and Vb as ~ 5 V. The c p value 
depends on the length of the wires that carry the sig- 
nal out of the cryostat, typically it is of order 400 pF. 
The amplifier gain, the Bessel frequency bandwidth, the 
duration of the acquisition window and the sampling fre- 
quency are set typically at 5000 V/V, 12 Hz, 5.008 s and 
125 Hz respectively. 

A typical signal acquired by the ADC, is shown in 
Fig. The signal was generated by a 2615 keV 7 particle 
fully absorbed in the crystal. The baseline voltage of the 
pulse is related to the thermistor temperature before the 
particle interaction, and the amplitude is related to the 
energy released. 
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Figure 3. Pulse generated by a 2615 keV 7 particle. The 
baseline is related to the temperature of the thermistor before 
the particle interaction, the amplitude carries information on 
the amount of energy released. 



III. NON-LINEARITIES OF THE RESPONSE 
FUNCTION 

The data analyzed in this paper have been obtained by 
exposing the bolometer to a 232 Th source. This source, 
together with an a line generated by 210 Po contamination 
of the crystal, allows the analysis of an energy range up 
to 5407 keV (see the energy spectrum in Fig. 0]). 
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Figure 4. Energy Spectrum. The a line at 5407 keV is gen- 
erated by 210 Po contamination of the crystal, the rest of the 
spectrum is due to 7's and /3's generated by a 232 Th source. 
The most prominent peaks are labeled with their energy (in 
keV). 



In this range several non-linearities of the response 
function are observed: 

1. The shape of the pulse is not constant with energy. 
The rise and the decay time, computed as the time 
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difference between the 10% and the 90% of the trail- 
ing edge and the 90% and 30% of the leading edge 
respectively, are correlated with energy (see Fig. EJ. 
This effect makes any pulse shape discrimination in 
the data analysis difficult. 
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Figure 5. Pulse shape parameters versus energy. The rise time 
(a) and the decay time (b) are anti-correlated and correlated 
with energy respectively. 

2. The amplitude of the pulse depends on the base 
temperature. In figure|5]the amplitude of 5407 keV 
events is plotted with respect to the pulse baseline, 
showing a clear anti-correlation. The temperature 
variation of the baseline is about 40 /j,K, which cor- 
responds to an energy variation of the crystal of 
about 400 keV. The pulse amplitude changes by 
about 1%, worsening the resolution of the detector. 
This effect can be partially corrected by applying 
a calibrated heat source to the detecto r 12 ' 13 and 
stabilizing the base temperature 14 . 

3. The amplitude dependence on energy is not lin- 
ear, and is usually parameterized with a polynomial 
function. The deviation from linearity of the data 
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Figure 6. Amplitude of 5407 keV pulses versus baseline. A 
change in the bolometer temperature also changes its re- 
sponse, worsening the resolution. 



is estimated using a linear calibration function 

Energy = constant • Amplitude (2) 

and then the difference between the energy esti- 
mated with this function and the nominal one is 
computed. The residuals evaluated on the peaks 
labeled in Fig.g]are shown in Fig. [7] The 5407 keV 
line is not included, since a particles have a quench- 
ing factor different from 7 and j3 particles^. If the 
true calibration function were a line, the residuals 
would be much smaller than the energy resolution, 
that is of order 5 keV FWHM. The data, however, 
clearly indicate that a line is not sufficient to de- 
scribe the calibration function. 



IV. THERMISTOR AND BIASING CIRCUIT MODELS 

The detector contains many elements that can generate 
the observed non-linearities. The heat capacitance and 
conductance of the crystal, of the Teflon supports and of 
the thermistor, all depend on the temperature. The crys- 
tal capacitance for example is not linear in temperature, 
as it follows the Debye law (C oc T 3 ). The problem of 
modeling the response function is that not all elements 
can be measured with satisfactory precision. In partic- 
ular the contact heat conductance between the crystal 
and its supports changes with the detector configuration 
and the heat conductance of the glue varies from bolome- 
ter to bolometer because of the weak reproducibility of 
the glue deposition. In principle the knowledge of all 
parameters is needed if we want to predict the response 
function. However, we can narrow our focus to the as- 
pects of the physical system which produce the dominant 
contributions to the non-linearity. In this section we will 
develop a model of the thermistor response and of the 
biasing circuit. It will be shown that they can account 
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Figure 7. Residuals obtained using a linear calibration func- 
tion. The difference AE between the estimated energy and 
the nominal one is not compatible with zero, considering that 
the energy resolution is ~ 5 keV FWHM. The error bars refer 
to the uncertainty on the estimated peak position, that de- 
pends on the FWHM and on the number of events N in the 
peak as FWHM/(2.35V^V). 



for all the observable non-linearities, without using any 
thermal model. 

First we have to extract the temperature of the ther- 
mistor from the measured resistance. Unfortunately we 
cannot use Eq. ([IJ, since parameters Rq and To cannot 
be measured with precision. This is due to the fact that 
when the thermistor is attached to the crystal, the glue 
stresses the device modifying its features in an unpre- 
dictable way. However, since we are interested in mea- 
suring energy depositions and not absolute temperatures, 
we can just evaluate the temperature variation. The re- 
sistance variation, AR, in terms of the temperature vari- 
ation, AT, can be approximated as: 



AR(AT) 



R(T + AT) - R(T) 
R(T) [exp(-77AT/T) 



where: 



n 



dlogR 



dlogT 



= 7I0. 



„ gen 

Rq 



(3) 



(4) 



is the sensitivity of the thermistor and its value is of order 
10. The reason for expanding the exponential argument 
in terms of AT/T instead of the entire formula is that 
it provides a better approximation. The approximation 
is valid for AT/T -C 1, or, in terms of particle energy, 
for E <C CT. Since CT is of order 100 MeV (C and 
T are about lOMeV/mK and 10 mK respectively), the 
approximation holds in the entire energy ranged . The 
advantage of Eq. ([3]) is that it allows AT to be extracted 
from AR using a measurable parameter (the resistance 
i?(T)) and that the unmeasurable parameter (rj/T) is 



just a scale factor. In order to simplify the notation we 
define a = rj/T for future calculations. 

In the next section we will put in evidence the non- 
linearities generated by Eq. ([3]) making the following as- 
sumptions: 

1 . The temperature response of the detector after an 
energy deposition E can be parametrized as 



AT(t,E) = A{E)U{t) 



(5) 



2. 



where U(t) is a function representing the time evo- 
lution of the temperature pulse and A{E) is the 
amplitude of the pulse, which depends on the en- 
ergy deposited. 

The pulse shape U(t) is independent of the energy 
and has unitary amplitude. 

3. The amplitude A(E) is proportional to the energy. 

4. Both A and U are independent of the baseline tem- 
perature (the temperature before the energy depo- 
sition). 

The properties of the resistance signal, AR(t), obtained 
under these assumptions will be compared with the mea- 
sured voltage signal, assuming that the biasing circuit 
does not add non-linearities. 

The dependence of the shape parameters on energy 
(see Fig. [SJ is obtained computing the time derivative of 
the signal AR(t,E) normalized to its amplitude Ar(E); 

AR'(t,E) _ exp[-aA(E)U(t)]aA(E)U'(t) 



Ar(E) 



U'(t) [1 



exp[-aA(E)} 
- aA(E) U(t)} 



(6) 



where Ar{E) — mm t [AR(t, E)] (since AR(t, E) is a neg- 
ative signal) and U'(t) is the time derivative of U (t). This 
expression indicates that the signal is slower when the en- 
ergy increases, the same effect seen on the decay time of 
the pulse. This equation does not explain the behavior 
of the rise time, that is dominated by the biasing circuit 
(see later in the text). 

The dependence of the pulse amplitude on baseline is 
also straightforward to derive. The resistance variation of 
the pulse baseline depends on the temperature according 
to the expression: 



AR B = R(T) [exp(-aAT s ) - 1] 



(7) 



where AT^ is the temperature variation with respect to 
T. When energy is released, the temperature increases 
by A(E), and the total resistance change will be: 

Ai? Tot = R(T) [exp(-a(AT B + A(E))) - 1] . (8) 

The amplitude of the pulse, however, is computed with 
respect to the baseline. In terms of resistance variation 
it will be: 



A R (E, AR b ) = AR Tot - AR B 
AR b 



1 



R(T) 



Ar(E,AR b =0). (9) 
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Recalling that a baseline increase in Volts corresponds to 
a resistance decrease, this expression describes the anti- 
correlation ol the pulse amplitude with its baseline (see 
Fig. HI). 

The non-linearity of the calibration function comes di- 
rectly from Eq. ([3]). If we assume that A(E) = E/C, 
the pulse amplitude Ar(E) at the second order increases 
with energy as: 



A R (E) = R(T) 



E 

- a c 



E 



(10) 



Recalling that the signal in Volts has opposite sign of the 
resistance variation, this is the same behavior of Fig. [7J 
The thermistor formula ([3]) explains at this point at 
least qualitatively all the data, except for the rise time of 
the pulse. The calculations showed in fact that the shape 
parameters should get slower with amplitude. The rise 
time, however, becomes faster (see Fig. [SJ. This effect 
can be generated by the biasing circuit, whose parasitic 
capacitance, c p , introduces a cutoff that varies with the 
pulse amplitude. We solve the circuit of Fig. [2] approxi- 
mately, using Rl >• R and Rl ~> AR(t). We have that 
the voltage signal AV(t) in terms of the resistance signal 
AR(t) can be expressed as: 

AV(t) + [R(T) + AR(t)} c p AV'(t) = -^-AR(t), (11) 

where AV(t) = V{t) — Vr, and Vr is the voltage 
corresponding to the static resistance R{T): Vr — 
-V B R(T)/[R(T) + R L }. Expression (11]) has the form of 
a low pass filter whose cutoff, [c p (R(T) + Ai?(t))] _1 , de- 
pends on AR(t) and hence on its amplitude. The higher 
is the absolute value of the amplitude, the higher is the 
cutoff (since the amplitude is negative) , making the volt- 
age pulse faster. Since R(T)c p is of order 30 ms, this 
effect is visible only on the fast part of the voltage pulse 
(the rise time) and competes with the thermistor slowing 
down. 



V. LINEARIZATION OF THE RESPONSE FUNCTION 



where AV(t) is defined in the previous section. Then 
from Eq. the quantity AS(t) = a AT(t) is computed: 



AS(t) = -log 



AR{t) 
R{T) 



(13) 



The advantage of this algorithm is that it depends on 
parameters that can be measured with sufficient preci- 
sion. These parameters are the static bolometer resis- 
tance R(T) and the voltage across it Vr, the load re- 
sistance Rl, the parasitic capacitance c p and the bias 
voltage Vb- The unknown quantities have been hidden 
in a common scale factor (cr), that is not relevant if we 
are only interested in removing the non-linearities. This 
algorithm has been implemented as a digital filter on the 
voltage data sampled by the ADC. The derivative of AV 
has been computed as AV- = (V i+ i— V;_i)/(2 At), where 
i is the sample index and At is the sampling period. A 
comparison between the original AV(t) signal and the 
same signal transformed into AR(t) and AS(t) is shown 
in Fig. E 
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It has been shown that the behavior of the data can 
be explained assuming that the thermal response AT(t) 
does not depend on the energy released or on the vari- 
ation of the baseline bolometer temperature. The non- 
linearities of the voltage signal seem to be generated by 
the thermistor and by the biasing circuit. Here we derive 
an algorithm to extract a signal proportional to AT(t), 
and we will test its performance on data. 

First the variation of the thermistor resistance is ex- 
tracted from V(t) solving the biasing circuit of Fig. [2] 

= AV(t)[R(T) + R L ] + R L R(T) c p AV'(t) 

[ ' AV{t) + R L c p AV'{t) + V B Rl/[R{T) + R L ] ' 

(12) 



Figure 8. Comparison between the original AV(t) signal gen- 
erated by a 2615 keV 7 particle, and the same signal trans- 
formed into AR(t) and AS(t). The AR(t) signal is faster 
than the AV(t) one, because the low pass filter of the biasing 
circuit has been deconvoluted. The AS(t) signal is even faster 
since the thermistor slowing down effect has been removed. 

The AS" data have been analyzed with the same anal- 
ysis procedure applied on AV data. The features of the 
transformed signal are: 

1. The shape of the pulse has a much weaker depen- 
dence on energy. The correction is particularly ev- 
ident on the decay time, where the correlation has 
been reduced by a factor 6 (see Fig. 
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Figure 9. Pulse shape parameters versus energy. The AS data 
(light gray) are less dependent on energy than AV data (dark 
gray), the correction is particularly effective on the decay time 
(b). The rise time (a) of AS data is faster, because the low 
pass filter of the biasing circuit has been deconvoluted. 



2. The amplitude dependence on the baseline is re- 
moved. The slope is reduced by a factor 20 (see 

Fig. nun. 
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Figure 10. Amplitude of 5407 keV pulses versus baseline. The 
AS data (light gray) are less correlated with the baseline than 
AV data (dark gray). The ratio of the slopes is about 20. 
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Figure 11. Residuals of the linear calibration function. The 
residuals of AS data (light gray) are very close to zero and 
have no trend compared to AV data (dark gray). 
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In conclusion, the model we have developed describes 
nearly all the observable non-linearities of the response 
function. It relies only on fixed parameters, that can be 
measured on the detector. The data transformed with an 
algorithm derived from the model provide a much more 
linear response function, improving the operation and the 
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